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ABSTRACT 

We report the detection of HCO + ( J=4— >-3) emission in the Cloverleaf Quasar at z=2.56, using 
the IRAM Plateau de Bure Interferometer. HCO + emission is a star formation indicator similar 
to HCN, tracing dense molecular hydrogen gas (rt(H2) ~ 10 5 cm~ 3 ) within star-forming molecular 
clouds. We derive a lensing-corrected HCO + ( J=4— »3) line luminosity of ^hco+(4-3) = ^ ^-3) x 

10 9 (/iL/H) _1 K kms _1 pc 2 , which corresponds to only 48% of the HCO + ( J=l— >0) luminosity, and 
<4% of the CO( J=3— >2) luminosity. The HCO + excitation thus is clearly subthermal in the J=4— s>3 
transition. Modeling of the HCO + line radiative transfer suggests that the HCO + emission emerges 
from a region with physical properties comparable to that exhibiting the CO line emission, but 2x 
higher gas density. This suggests that both HCO + and CO lines trace the warm, dense molecular 
gas where star formation actively takes place. The HCO + lines have only ^2/3 the width of the 
CO lines, which may suggest that the densest gas is more spatially concentrated. In contrast to the 
z=3.91 quasar APM 08279+5255, the dense gas excitation in the Cloverleaf is consistent with being 
purely collisional, rather than being enhanced by radiative processes. Thus, the physical properties of 
the dense gas component in the Cloverleaf are consistent with those in the nuclei of nearby starburst 
galaxies. This suggests that the conditions in the dense, star-forming gas in active galactic nucleus- 
starburst systems at early cosmic times like the Cloverleaf are primarily affected by the starburst 
itself, rather than the central active black hole. 

Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Investigating the dense molecular interstellar medium 
ON i (ISM) in distant galaxies is of fundamental importance 
ON ■ for our general picture of galaxy formation and evolution 
' in the early universe, as it is found in the regions where 
active star formation occurs. Due to the fact that CO 
exhibits the brightest emission lines of all molecules, it 
is a good tracer for molecular clouds and the diffuse, 
j— ( . gaseous ISM; i.e., the total amount of potent ial fu el for 
I* ■ star formation (see Solomon & Vanden Bout 120051 for a 
.£h ' review). However, the low density required to excite CO 
^ . (> 300 cm -3 ) also means that it is not a specific tracer of 
the molecular cloud cores where star formation actively 
C$ ' takes place. In contrast, recent studies have shown that 
high dipole moment molecules like HCO + and HCN are 
substantially better tracers of cloud cores. This is due 
to the fact that such molecules trace much denser gas 
(> 10 5 cm~ 3 ) than CO emission (e.g., Gao & Solomon 
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I2004ail2004bh . 

As they only trace the densest regions of the ISM (and 
are less abundant), emission from rotational transitions 
of dense gas tracers is typically by at least an order of 
magnitude fainter than emission from CO lines. Thus, 
studies of the dense ISM in high redshift galaxies are cur- 
rently focused on only few, bright targets. Consequently, 
high-z HCN and HCO + emission was hrst detected in the 
Cloverleaf quasar (2=2.56), the brightest known CO line 
em itter at the time (Solomon et al. 120031 : Riechers et 
al. I2006f ). Both molecules were detected in the ground- 
state J=1^»0 transitions, which trace the full amount of 
dense gas. Despite a number of physical and chemical 
processes that affect HCN and HCO + in different ways, 
their respective J=l— >0 lines show comparable strengths, 
indicating that the emission is likely optically thick. 

To understand the physical conditions in the dense, 
star-forming molecular gas in more detail, it is nec- 
essary to study multiple rotational lines from dense 
gas tracers. Previous searches for HCO + ( J=4— s>3) and 
HCN(J=4— >3) line emission in the C loverleaf quasar 
were unsuccessful (Wilner et al. 119951 : Solomon et al. 
120031: see also Barvainis et al. 119971 ). The limits indi- 
cate a comparatively low excitation of HCN, but are not 
constraining for HCO + , given the measured J=l— >0 line 
luminosities. 

In this Letter, we report the detection of 
HCO + (J=4— >3) emission toward the Cloverleaf quasar 
(2=2.56), using the IRAM Plateau de Bure Interferom- 
eter (PdBI). This enables us to investigate the dense 
molecular gas excitation in a high-z galaxy, constraining 
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Fig. 1. — Velocity-integrated PdBI map of HCO+(J=4— >3) line emission over 284 kms -1 (approximately the line FWHM) toward the 
Cloverleaf. At a resolution of 5.3"x4.7" (as indicated in the bottom left), the emission remains unresolved (maxi mum lens image separation: 
~1.4"). The cross indicates the center of the CO(J=7— >6) emission (color scale in the right panel; Alloin et al. [19971 ). Contours are shown 
in steps of lc=0.31 mjy beam - 1 , starting at ±2<r. 

the physical conditions for star formation out to the 
early universe. We use a concordance, fiat ACDM 
cosmology throughout, with Hp= 71 kms" 1 M pc -1 , 
fi M =0.27, and S1 A =0.73 (Spergel et al. 12001 [2001 . 



2. OBSERVATIONS 

We observed the HCO + (J=4^»3) transition line 
(V re8t = 356.734288 GHz) toward H1413+117 (the 
Cloverleaf quasar), using the PdBI. A t the target red- 
shift of z=2. 55784 (e.g., Barvainis et al. 119941: Weifi etal. 
I2003D . the line is shifted to 100.267GHz (2.99mm). A 
total bandwidth of 580 MHz (^1700 kms -1 ) was used to 
cover the HCO + (J=4^»3) line and the underlying 3 mm 
(rest-frame 840 /im) continuum emission. Observations 
were carried out under acceptable to good 3 mm weather 
conditions during 10 tracks in D configuration between 
2006 May 12 and August 20. The total integration time 
amounts to 70.2 hr using 4, 5, or 6 antennas, result- 
ing in 10.8 hr 6 antenna-equivalent on-source time after 
discarding unusable visibility data. The nearby sources 
B1354+195 and B1502+106 (distance to the Cloverleaf: 
9.0° and 12.0°) were observed every 20 minutes for point- 
ing, secondary amplitude and phase calibrations. For 
primary flux and bandpass calibration, several nearby 
calibrators (MWC349, CRL618, 3C84, 3C273, 3C279, 
3C345, and 3C454.3) were observed during all runs, lead- 
ing to a calibration that is accurate within 10-15%. 

For data reduction and analysis, the IRAM GILDAS 
package was used. All data were mapped using the 
CLEAN algorithm and 'natural' weighting; this results in 
a synthesized beam of 5.3 "x 4.7" (~40kpc at z = 2.56). 
The final rms is 0.31 mJy beam -1 over 95 MHz (cor- 
responding to 284 kms -1 ), and 0.79mJy beam -1 over 
15 MHz (45 kms -1 ). 

3. RESULTS 

We have detected HCO + (J=4— >3) line emission to- 
ward the Cloverleaf quasar (z—2.56) at 6er signifi- 
cance (Fig. [IJ. From Gaussian fitting to the line pro- 
file, we obtain a HCO + (J=4->3) line peak strength 



of 1.77±0.35mJy at a FWHM of 288±79 kms -1 , on 
top of 0.29±0.14mJy continuum emission (consistent 
with 0.3±0.1mJy as measured at 93 GHz an d the 
dust spectra l ener gy distribution; Henkel et al. 120101 : 
Weifi et al. [2003). This corresponds to a velocity- 
integrated emission line strength of 0.54±0.09 Jy kms -1 , 
and a line luminosity of £hco+(4-3) = (1-6 ± 0.3) x 
10 9 (/xl/11) -1 K kms -1 pc 2 (where \i\=\\ is the le nsing 
magnification factor; Venturini & Solomon [2J303J) , i.e., 
only <4% of the CO(J=3->2) luminosit y, and ~35% of 
the C N(iV= 3->2) luminosity (Weifi et al. 120031: Riechers 
et al. l2007al) . 

The line FWHM is only ~7 0% o f that of the 
CO(J=3^2) line (Weifi et al. [2003). With this 
new constraint, we re-visit ed th e HCO + (J=l— >0) line 
data by Riechers et al. (2006; taken with substan- 
tially narrower bandwidth), which already showed 
some evidence for a narrower line than CO. The 
HCO + (J=1^0) data can be fitted well with a 
(beam-corrected) line peak strength of 0.249±0.028mJy 
at a FWHM of 263±52 kms -1 This corre- 
sponds to a velocity-integrated emission line strength 
of 0.069±0.008 Jy kms -1 , and a line luminosity of 
L hco+(i-o) = (3-3 ±0.3) x 10 9 (Ai L /ll) _1 K kms -1 pc 2 , 
which is consist ent w ith the previously derived value 
(Riechers et al. 2006) within the errors. More im- 
portantly, this corresponds to a HCO + J=4— >3/1^0 
line brightness temperature ratio of r4i=0.48±0.11, i.e., 
the HCO + ( J=4— s>3) line is clearly subthermally ex- 
cited0 This is cons istent with the previ ous li mit of 
r 4i <4 (Wilner et al. [19951: Riechers et al. (2006). We 
also set a 3a lower limit of r 4 (HCO+/HCN)>0.59 on 
the HCO+/HCN J=4->3 line ratio, consistent with the 
J =l— >Q line ratio of ~0.8 within the errors (Rie chers e t 
al. 120061: HCN J=4^3 limit from Solomon et al. [200l . 

8 This width is also consistent with that of the HCN(J=1— >0) 
line (Solomon et al. l2003l) . 

9 Or optically thin, but see Sect. 4. 
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Fig. 2.— Left: PdBI HCO+( J=4->3) spectrum of the Cloverleaf 
at 15 MHz (45 kms -1 ) resolution (histogram), along with a Gaus- 
sian fit to the line and continuum emission (black curve). The line 
FWHM is 288±79 kms" 1 . The velocity scale is relativ e to the 
source 's redshift of 2=2.55784 (e.g., Barvainis et al. [19941 ; Weifi et 

ai. rsooiy 

4. HCO+ AND CO LINE EXCITATION MODELING 

Based on the HCO + excitation ladder of the Clover- 
leaf, we can constrain the line radiative transfer of the 
dense molecular gas component through Large Velocity 
Gradient (LVG) models, treating the gas kinetic temper- 
ature and density as free parameters. To maximize the 
available observational constraints on these parameters, 
we here simultaneously model the HCO + and CO exci- 
tation, and also require that our solutions are consistent 
with the dust sp ectral energy distribution of the Clover- 
leaf (Weifi et al. [2003|). 

Our models use the H CO + and CO collision rates 
from Flower (|1999t 120011 ). We adopt a HCO+ abun- 
dance per velocity gradient of [HCO + ]/(dt>/dr) = 1 x 
lO^pcQtms- 1 )- 1 , and [HCO+/CO]=10- 4 (e.g., Wang 
et al. 120041 ). The HCO+ data are fit well by a spher- 
ical, single-component model with a kinetic tempera- 
ture of Tkm=50K, a gas density of n gas =10 4 8 cm~ 3 , 
and a CO disk filling factor of 22% (rco=785pc, as- 
suming /il=11; right panel in Fig. [3]). In this model, 
the HCO + J=l— >0 and 4h>3 lines have optical depths 
of Ti_o=3-l and T4_3=25.4 and excitation temperatures 
of T4" =36.7K and T C 4 - 3 =24.4K. This suggests that 
the emission in both transitions is optically thick, and 
that the HCO + (J=4— >3) line, indeed, is subthermally 
excited. 

The (literature) CO data are fit well by a spheri- 
cal, single-component model with a kinetic temperature 



of Tkin=50K, a gas density of n g 



= 10 



4.5 , 



and a 



— 5.5x higher surface filling factor than HCO + (left panel 
in Fig. [3]). These parameters are comp atible with the CO 
excitation analysis of Bradford et al. (|2009T ) . Our anal- 
ysis thus suggests that the CO and HCO + excitation in 
the Cloverleaf can be modeled simultaneously, with the 
same kinetic temperatures and only a factor of 2 differ- 



ence in gas densities. 

To explore the remaining uncertainties, we also at- 
tempted to fit the HCO + data with Tki n and n gas fixed to 
those of the CO model, but varying the relative molecular 
abundance. We find an acceptable solution when increas- 
ing the relative HCO + abundance to [HCO + ]/(di;/dr) = 
1 x 10- 8 5 pc(kms~ 1 )- 1 , i.e., [HCO+/CO] = 10- 3 - 5 , and 
the HCO + surface filling factor to 24% (see dashed and 
dotted lines in Fig. [3]). 

Overall, the models suggest a relatively high median 
gas density in this galaxy, and that the CO and the 
HCO + emission likely trace the same warm, dense molec- 
ular ISM phase, with HCO + tracing the densest ~15- 
20% of the gas. Given that the HCO + lines have only 
~2/3 of the width of the CO lines, this may suggest that 
the emission region with the densest gas is more spatially 
concentrated than the overall CO emission, such as, e.g., 
in a nuclear starburst. Higher spatial resolution HCO + 
observations are required to confirm this scenario, and to 
investigate potential differential lensing effects. 

5. DISCUSSION 

We have detected bright, but subthermally excited 
HCO + (J=4^3) emission toward the Cloverleaf quasar 
at z=2.56. Based on excitation modeling, we find that 
the warm, dense gas traced by HCO + appears to be asso- 
ciated with the warm gas phase traced by the CO lines, 
picking out its densest regions. 

The HCO + excitation in the Cloverleaf is consistent 
with that seen in the starburst nucleus of NGC 253, with 
HCO+ J=4->3/ l->0 r atios of r 4 i=0.48±0.11 and 0.53 
(Knudsen et al. 120071 ). respectively. Intriguingly, the 
Cloverleaf and the nucleus of NGC 253 also have com- 
parable HCO+/HCN J=1^0 line ratios of -0.8. The 
m in the Cloverleaf is higher than that in the infrared- 
luminous galaxies NGC 6240 (0.21±0.06) and Arp220 
(0.33±0.12; Greve et al. I2009D . However, this is likely 
due to the fact that the line ratios are averaged over vir- 
tually the entire molecular line emission regi ons, r ather 
than just the nuclei. As shown by Iono et al. (|2007t their 
Fig. 10), the T4i in NGC 6240 scatters up to values of 
~0.6 within the HCO + -emitting region, indicating that 
the ratio is comparable to the Cloverleaf in the densest 
regions. Thus, it seems plausible that the dense gas ex- 
citation in the Cloverleaf is comparable to what is found 
in the nuclear regions of nearby starburst galaxies and 
luminous infrared galaxies. 

The relatively high median gas density in the Clover- 
leaf suggested by the HCO + and CO observations is also 
consistent with its location on th e HCO + -far-infrared lu- 
minosity relation (Riechers et al. 2006) within the frame- 
work ofthe model description of Krumholz & Thomp- 
son (J2007) . In fact, it may be the most direct evidence 
that the increasing slope in dense gas-star formation rela- 
tions observed towa rd the most luminous, high r edshift 
systems (Gao et al. I2007t Riechers et al. I2007b0 is in- 
deed related to an elevated median gas density relative 
to lower-luminosity systems. 

Besides APM 08279+5255 (2=3.91), the Cloverleaf is 
only the second high-z galaxy in which multiple tran- 
sition s of a dense gas tracer we re det ected (e.g., Wagg 
et al. [20051 Garcia -Burillo et al. [20061: W eifi et al. l2007t 
Guelin et al. I2007t Riechers et al. 120091 12010D . Model- 
ing of the HCN line ladder in APM 08279+5255 suggests 
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Fig. 3. — C O (le ft) and HCO+ (right) exci tation ladders (points) and LV G mod els (lines) for the Cloverleaf. Left: The CO data (circle: 
Tsuboi ct al. 1999;; triangles: Barvainis et al. 119971 : filled square: Weifi et al. 120031 : filled circles: Bradford et al. 20Q3) are fit well by a gas 
component with a kinetic temperature of T^ ln =50K and a gas density of p gas =10 4 - 5 cm" 3 . Right: The HCO+ data are fit well by a dense 
gas component with Tj i j n =50K and p gas =10 4 ' 8 cm -3 (i.e., 2xp gaa of CO; solid line). For comparison, the dashed line shows a model with 
the same p gas as in the CO model. The dotted line shows a model with the same p gas as in the CO model, and a 10°- 5 x higher HCO+ 
abundance as in the other models. All models are scaled to the same HCO + (J=l— >0) flux. 



that the emission in high- J HCN transitions is substan- 
tially enhanced by radiative excitation through pump- 
ing of mid-infr ared r o- vibrational lines (Weifi et al. I2007t 
Riechers et al. 120101 ) . In contrast, the HCO + excitation 
in the Cloverleaf is consistent with purely collisional ex- 
citation. Given the comparable critical densities of HCN 
and HCO + , this suggests that we have identified a clear 
difference in the dense gas excitation conditions between 
these two high-z systems. 

This investigation highlights the importance of study- 
ing the excitation of dense gas tracers to understand 
differences in the conditions for star formation in high 
redshift galaxies. Such studies will become routine with 
the advent of broad instantaneous bandwidth systems as 



part of future facilities such as the Atacama Large (sub-) 
Millimeter Array (ALMA) , which will allow to frequently 
cover lines of multiple dense gas tracers as part of 'stan- 
dard' high-z CO observations. 
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